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■ ABSTRACT 

o • 

\^ BL Lacertae has been the target of several observing campaigns by the Whole Earth Blazar Telescope (WEBT) collaboration and is one of the 
best studied blazars at all accessible wavelengths. A recent analysis of the optical and radio variability indicates that part of the radio variability 
is correlated with the optical light curve. Here we present an analysis of a huge VLBI data set including 108 images at 15, 22, and 43 GHz 
obtained between 1995 and 2004. The aim of this study is to identify the different components contributing to the single-dish radio light curves. 
' We obtain separate radio light curves for the VLBI core and jet and show that the radio spectral index of single-dish observations can be used 
O ^ to trace the core variability. Cross-correlation of the radio spectral index with the optical light curve indicates that the optical variations lead the 
I radio by about 100 days at 15 GHz. By fitting the radio time lags vs. frequency, we find that the power law is steeper than expected for a freely 
O , expanding conical jet in equipartition with energy density decreasing as the square of the distance down the jet as in the Konigl model. The 
• analysis of the historical data back to 1968 reveals that during a time range of 16 years the optical variability was reduced and its correlation 
^ with the radio emission was suppressed. There is a section of the compact radio jet where the emission is weak such that flares propagating 
. • down the jet are bright first in the core region with a secondary increase in flux about 1.0 mas from the core. This illustrates the importance 
^ ' of direct imaging to the interpretation of multi-wavelength light curves that can be affected by several distinct components at any given time. 
>v>( ^ We discuss how the complex behaviour of the light curves and correlations can be understood within the framework of a precessing helical jet 
, model. 

Key words, galaxies: active - galaxies: BL Lacertae objects: general ~ galaxies: BL Lacertae objects: individual: BL Lacertae - galaxies: jets 
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1. Introduction 

The term "blazars" identifies a family of radio-loud active 
galactic nuclei (AGN) showing strong variability at all wave- 
lengths from radio to y-rays, high degrees of polarization, and 
apparent brightness temperatures exceeding the Compton limit 
(e.g., Urrv 1999). During the last decades a rather general con- 
sensus on the global mechanism responsible for the emission 
has been achieved: a supermassive black hole surrounded by a 
massive accretion disk feeding a powerful jet closely aligned 
to the Une of sight. Relativistic electrons in the jet plasma pro- 



Send offprint requests to: U. Bach, e-mail: bach@to . astro . it 

* For questions regarding the availability of WEBT data please con- 
tact Massimo Villata (e-mail: villataOto . astro . it) 



duce the soft synchrotron photons (from radio wavelengths to 
the UV band and sometimes even X-rays) through, while hard 
photons (from X-rays to y-rays) are like ly produced by inverse 
Compton scattering (e.g. JSambruna et~al..,1996t and references 
therein). Blazars are divided into two subclasses: flat-spectrum 
radio quasars and BL Lac objects, whose main difference is the 
lack or weakness of emission lines in BL La c objects (strongest 
rest -frame emission line w idth < 5A; e.g.. lStickel et alJll99ll 
and lUrrv & Padovanil 19951 for a review). 

BL Lacertae, located in the nucleus of a moderately brigh t 
elliptical galaxy at a redshift of z = 0.069 (Miller et alJll978h 
is the prototype of the class of BL Lac objects. However, emis- 
sion lines with equivale nt width in excess of 5 A have been 
detected several times dVermeulen et alJ Il995t ICorbett et alJ 
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Il996l I2OOOI) . suggesting that BL Lac objects can also have 
a broad line region (BLR), but it would be usually out- 
shone by the beamed synchrotron emission of the jet. BL Lac 
has been studied intensively since its discovery and was the 
target of several mult i-wavelengths campaigns (Bloom et aL 
1997 |ISambruna et al Il 999: Madeisk i etalJl999^lviUata et all 
2002^. 'Bottcher et al."2 003ii: iVillata et al.l l2004a*b^. In particu- 
lar, the data collected during the Whole Earth Blazar Telescope 
(WEBT)' campaigns provide an unprecedented time sam- 
pling for radio and optical light curves from 1994 up to 
now ( Villata et al. 2002, 2004a b). Cross-corelation analysis of 
these light curves revealed well correlated variability in the 
radio bands, where variations at higher frequencies lead the 
lower-frequency ones by several days to a few months depend- 
ing on the frequency separation. The detection of a fair corre- 
lation between the optical R band variations and the ratio be- 
tween the 22 GHz and 5 GHz radio flux densities suggests that 
the optical emission and part of the radio emission arises from 
a common origin in the inner portion of the jet. 

Very Long Baseline Interferometry (VLBI) observations of 
the parsec-scale structure of BL Lac reveal a brig ht compact 
radio core and a jet extending to the south (e.g., Mutel et all 
1 1990) . The jet emerges at a position angle of P.A. x 195° and 
displays a small bend at ~ 4 mas distance from the core towards 
P.A. ^ 160°. A number of superluminal jet components have 
been obser ved displayirig bent trajectories and speed s from 3 c 
to 9 c [e.g.. iMutel et al.lll990t JPenn et alJ:2 000: St irUng et alJ 



l2003t iKellermann et alJT2004t |jorstad'etalJl2005<) . Evidence 
for precession of the jet nozzle arising from periodic varia- 
tions of the electric vector position angle (EVPA) at millimetre 
wavelengths and the position angle of newly emerging jet fea- 
tures ( Stirling et al. 2003) are still under debate (Mutel & DenJ 
I2OO5I) . A comparison of the position angles of optical polariza- 
tion vectors with the orientation of the polarization vectors in 
5 GHz VLBI images suggests that the optica l emission is more 
likely to arise from the core than from the jet JGabuzda & Sitkol 

Here we present an investigation of the flux density evo- 
lution of the parsec- and subparsec-scale structure of BL Lac 
at 15 GHz, 22 GHz, and 43 GHz using a dataset of 108 VLBI 
observations performed between 1995 and 2004. The aim of 
this study is to disentangle the dififerent variability patterns that 
seem to be present in the single-dish radio light curves and 
to identify the region of the radio events that seem to corre- 
late with the optical variability. Throughout this paper we will 
assume a flat universe model with the following parameters: 
Hubble constant Hq = 71 kms"' Mpc"', a pressure-less mat- 
ter content - .3, and a cosmological constant Q^i - 0.7 
iSnergel et al. P2003'). At the redshift of BL Lac this corresponds 
to a luminosity distance of 308 Mpc and to an angular resolu- 
tion of 1.3pc/mas. Cosmology-dependent values quoted from 
other authors are scaled to these parameters. 

The observations and data reduction procedures will be de- 
scribed in Sect. 121 Our results will be presented in Sect. |3l and 
discussed in Sect. |2 At the end we will give a summary with 
the conclusions in Sect.|5l 



2. Observations and data reduction 

2.1. Radio and optical light curves 

The optical and single-dish radio flux densities used in this 
work have already been presented and analysed bv lVillata et alJ 
(2004a b). The radio data at 4.8 GHz, 8.0GHz, and 14.5 GHz 
(labelled as 5, 8, and 15 GHz throughout the article) were ob- 
tained by the AGN monitoring at the Universit y of Michiga n 
Radio Astronomy Observatory (UMRAO'; AUe r et alJl2003l) . 
while the 22 GHz and 37 GHz flux densities were measured 
with the radio telescope of the Metsahovi Radio Observatory 
(Teriisrantaet al. 1998, 2004, 2005). The optical R band light 
curve was obtained by assembling data from various observa- 
tories of the WEBT collaboration; details on the light curve 
construction and th e subtraction of the host ga laxy contribu- 
tion can be found in Villata et alJll2002l ' 2004a b). Here we use 
these data for a comparison with the structural evolution of the 
VLBI core and jet structure on parsec and subparsec-scales. 



2.2. VLBI data 

Most of the VLBI data used in this study were provided by 
the authors from various large observing campaigns, where 
BLLac was observed as a target source^* or as a polarization 
calibrator'*. A list showing the beam size, total flux density, 
peak flux density, rms noise level, and the references where 
the data was first published for all our VLBI images is given 
in TablefTlravailable in the electronic version). For the details 
on the data reduction we refer to the references of the obser- 
vations. The data were provided as calibrated (u, v)-data files 
and no further self-calibration was applied. Three epochs (ref- 
erence 6 in Table Q are newly reduced observations where 
BL Lac was used as a polarization calibrator. Those were 
obtained with the full VLBA and the 100 m Effelsberg tele- 
scope at 15 GHz and were correlated at the VLBA correlator 
in Socorro, NM. The post-correlation analysis was done using 
NRAO's Astronomical Image Processing System (AIPS). After 
loading the data into AIPS, the standard amplitude and phase 
calibration was performed. Both imaging an d phase- and am - 
plitude self calibrati on was done i n Difmap dShepherdll 1 997l) . 
using the CLEAN (Hogb omI 19741) and SELFCAL procedures. 
All data were imaged with the same image parameters at each 
frequency using uniform weighting in Difmap. Images with 
matching cell-sizes and resolution (high frequency data were 
tapered to the lower frequency resolution) were produced for 
the spectral analysis. The analysis of the images, including 
spectral index maps, slices and flux density measurements, was 
done in AIPS. 

Sample images of BL Lac at 15 GHz, 22 GHz, and 43 GHz 
are shown in Fig.[n The jet is visible up to 10 mas (13 pc) from 



^ http : //www . astro . Isa . umich . edu/obs/radiotel/umr ao.html| 
^ VLBA 2 cm Surv ey y Cellermann et aL .1998; .Zensus et alj 



| 2002l: iKellerma nn et alj 12004)? MOJAVE jLister & Homan 



2005 ); 



http://www.to.astro.it/blazars/webt 



'Denn et al. ( 2000): 'Jo rstad et alJ <200U llOOSb ^ StMin^^^lTl2003l) : 

,Mutel&Denn (2005) 

Agudo et al. & Gomez et al. p riv. comm. target source published 

mlf 



a lGomez et alJj200dl200lLl2002l) 



U. Bach et al.: Structure and flux variability in tlie VLBI jet of BL Lacertae during the WEBT campaigns (1995-2004) 



3 



15 GHz 




22 GHz 




43 GHz 






- 




- 






















W 
















- 













-2 -4 



Right Ascension (mas) 



Fig.l. VLBI images of BL Lac. Left: 15 GHz image (epoch 2002.03) with a beam of 0.64 mas x 0.32 mas at P. A. -33°. Peak 
flux density is L4 Jy/beam and contours start at L5 mJy/beam, increasing by a factor of 2. Middle: 22 GHz image (2002.05) with 
a beam of 0.88 mas x 0.58 mas at PA. 27°. Peak flux density is L6 Jy/beam and contours start at 3mJy/beam, increasing by a 
factor of 2. Right: 43 GHz image (2002.05) with a beam of 0.58 mas x 0.34 mas at PA. 32°. Peak flux density is L3 Jy/beam 
and contours start at 3 mJy/beam, increasing by a factor of 2. 



the bright VLBI core and bends at about 4 mas from a position 
angle of ~ 195° to ~ 160° (measured counterclockwise from 
north). The VLBI core is clearly distinguished from the jet by a 
region of low emission (between 0.7 mas and 1.0 mas core dis- 
tance) most of the time. This becomes more evident in Fig.|2l 
where the source intensity profiles along PA. = 195° are shown 
(see next section for more details). 

3. Results 

3.1. VLBI light curves 

To quantify the variability of different regions of the source 
we extract light curves of the core and of diff'erent portions 
of the jet. The regions that we chose are displayed in Fig. |2l 
Based on the jet kinematics we expect that a flare which oc- 
curs in the core will propagate along the jet and cause the jet 
to brighten at a certain position as the feature passes by. We 
will later see that this is indeed the case. Therefore, the exact 
location of the flux density measurement along the jet seems 
not very important, since a flare that travels down the jet will at 
some time pass by. Remarkable in this context is the low emis- 
sion region around 0.7-1 mas from the core, where only very 
small variations occur even though the core and the jet around 
2 mas are highly variable and sev eral new jet components were 
reported during this tim e range (Denn et al. 2000; Jorstad^^^^ 
ioP.'S tirling et alJ200.3^ Kellermann et al..2004;.Jorstad et alJ 
2005). It is obvious that a measurement in this region would 
give completely diff'erent results than a measurement 0.5 mas 
farther down the jet. From Fig. |2] we can estimate the maxi- 
mum brightness that can be reached at a certain jet position. It 
seems that new components first rapidly fade as they separate 
from the core and then reappear at about 1 mas distance from 



the core. The maximum brightness is reached at about 1.5 mas 
and after that the components slowly fade as they travel down 
the jet. 

We have tested several alternative positions with different 
sizes to check how the light curves change at different core dis- 
tances. If we take measurements within a very small region, the 
integrated flux density goes down and the signal to noise ratio 
of the light curve decreases. On the other hand, if the region 
is too large, different jet components will overlap and smear 
out the variability. We find that the best choice is a small re- 
gion close to the core, where flares from the core might be de- 
tectable with only a short time delay, as well as two slightly 
larger regions after the low emission gap to trace the flux den- 
sity evolution farther downstream. 

The core is represented by a circular region of 0. 15 mas ra- 
dius around the image centre, the first jet region corresponds 
to a core separation of 0.2 mas < r < 0.7 mas, the second one 
to 1.5 mas < r < 2.5 mas, and the third one to 2.5 mas < r < 
3.5 mas. We have tested several methods to obtain the corre- 
sponding light curves, including 1) the task IMSTAT in AIPS, 
which integrates the flux density of the image within a pre- 
defined rectangular region, 2) model fitting of Gaussian com- 
ponents in Difmap; and 3) integrating the flux density along 
the slices presented in Fig |2] Among the alternatives, we find 
that a quite accurate and flexible method is to extract the light 
curves by summing the delta-function components within the 
specified regions from the final CLEAN models produced with 
Difmap. Another advantage of this method is that one can eas- 
ily change the limits of the regions and rerun the analysis. Since 
most of the data were reduced by various authors and the qual- 
ity of the data is inhomogeneous, we estimate the flux den- 
sity error to be 10%, which is typical for VLBI measurements. 
From our experience with the VLBA, 10% is a conservative 
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Fig. 3. Flux density evolution of different portions of the VLBI core-jet structure of BLLac at 15 GHz, 22 GHz, and 43 GHz 
compared with the 15 GHz single-dish light curve. Different symbols denote different separation from the core r (see text for 
more details). The grey box highlights an example of a nice case study, where there are three equally bright flares, but the VLBI 
data reveals that contrary to the others the middle one was dominated by the jet and not by the core. Note that due to the rapid 
fading of the jet we do not show a Jet3 curve at 43 GHz. 



assumption and the typical uncertainty at 15 GHz and 22 GHz 
might be as low as 5%. Our three VLBA+Eff'elsberg epochs 
(2002.03, 2002.51, and 2003.24) at 15 GHz have typical errors 
of 5% estimated from the gain corrections during the amplitude 
self-calibration. 

Figure |3] shows the final VLBI light curves in comparison 
with the 15 GHz single-dish light curve from UMRAO. One 
can clearly see that the single-dish light curve is usually domi- 
nated by variability of the VLBI core, but occasionally parts of 
the jet reach or even overcome the core brightn ess. The latter is 
most evident at 15 GHz before mid 1998. Villat a et all ( l2004al) 
calculated the flux density ratio (a "hardness ratio," which is 
comparable with the radio spectral index) between the 22 GHz 
and 5 GHz light curves and distinguished two dififerent kind of 
flares. Flares that are more pronounced at higher frequencies 
are called hard flares (inverted or flat radio spectrum), while 
flares that are dominant at lower frequencies are called "soft 
flares" (steep radio spectrum). 

Fro m a comparison of Fig. |3] with Fig. 8 in Ivillata et alJ 
(l2004al) . one can see that the hard flares are observed when 



the VLBI core is bright compared to the jet, while during soft 
phases the jet is relatively bright. For example, the second flare 
(1998.2) of the three equally bright flares in the single-dish 
light curve between 1997.7 and 1998.9 in Fig. El (grey box) 
is softe r than the adjacent ones (cf. also Fig. 8 iri mllata et all 
l2004ah . If we take a look at the VLBI light curves it becomes 
obvious that, although in the single-dish measurements the 
flares have a similar brightness, the 1998.2 flare has a large 
contribution from a brightening in the jet at about 2 mas core 
distance (jet2) and originates not only from the core. It seems 
therefore very important to know the jet contribution to the 
single-dish light curves if we intend to compare the radio vari- 
ations with X-ray or optical light curves, since they are most 
likely dominated by the emission from the jet-base and not 
from the outer optically thin jet emission. 

3.2. Radio spectral index 

The single-dish radio spectral index a (defined as Sy ^ v", 
where Sv is the flux density at frequency v) was calculated for 
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Fig. 2. Intensity profiles of the VLBI images of BLLac at 
22 GHz (top) and 43 GHz (bottom), combining all epochs. The 
figure illustrates the flux density behaviour of the jet ridge-line 
at a resolution of 0.5 mas (22 GHz) and 0.3 mas (43 GHz) at 
P.A. - 195°. Each curve represents one epoch between 1997 
and 2003. Vertical lines indicate the core and jet regions that 
were used to analyse the images. 
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all adjacent frequency pairs between 5 GHz and 37 GHz with 
separation shorter than 3 days. This limit is usually fulfilled for 
data from the same observatory (5, 8, 15 GHz from UMRAO 
and 22 and 37 GHz from Metsahovi) and therefore reduces the 
sampling significantly only for the 15 GHz to 22 GHz spectral 
index. The mean time sampling for the radio light curves is 
typically of 10 days, but after 1990 the time sampling becomes 
much better (~ 5 days). The spectral indices obtained for dif- 
ferent pairs of frequencies (5/8, 8/15, 15/22, and 22/37 GHz) 
all follow the same behaviour in time. In Fig. 0] (top panel) 
we show the spectral-index curve between 8 GHz and 15 GHz 
(cs/is), which is the best sampled one. One can see that most 
of the flares are accompanied by a flattening of the spectrum. 
The flattening seems to peak shortly before the flare reaches its 
maximum. A comparison between the 43 GHz VLBI core flux 
density and ag/is is shown in the bottom panels. The single- 
dish spectral index is able to trace the VLBI core variability in 
a very accurate way. We will test this using a discrete correla- 
tion function (DCF) in the next section. 

The time evolution of the spectral index of the VLBI core 
and the inner 3 mas of the jet can be seen in Fig. |5l Shown 
are the ridge-line profiles of BL Lac across the spectral index 
VLBI maps between 22 GHz and 43 GHz along P.A. 195°. The 
spectral index maps were produced by aligning the two images 



Fig. 4. Top panels: Radio light curve at 15 GHz in comparison 
with the spectral-index curve between 8 GHz and 15 GHz. One 
can see that most of the flares are accompanied by a flattening 
of the spectrum. The vertical dotted lines mark some examples 
and help to guide the eye. Bottom panels: Zoom into the time 
range of the VLBI monitoring at 43 GHz. The grey line repre- 
sents a three point average of the spectral-index curve. There 
is a close apparent correlation between the 43 GHz core flux 
density and o-g/ig (see also Fig.|6j. 



on the brightest component. To reduce the effect of a potential 
core shift between the two frequencies, we convolved the im- 
ages with a 0.7 mas circular beam. Since there are no strong 
spectral gradients visible in the core region that would indi- 
cate a core shift, we assume that the effect is negligible at the 
resolution of our images. The core spectrum in Fig.|5lis signif- 
icantly flatter than the jet for most of the time especially when 
it is flaring (e.g., epoch 1998.61). It therefore seems reasonable 
that the variability of the single-dish radio spectral index is a 
good tracer for the variations in the VLBI core. 
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Fig. 5. Ridge-line profiles of the spectral index maps between 
22 GHz and 43 GHz showing the core and the inner 3 mas of 
the jet along P.A. 195°. The original images were tapered to the 
same resolution and convolved with a circular beam of 0.7 mas. 



lights radio flares that are more pronounced at higher frequen- 
cies. On the other hand, the radio spectral index and the VLBI 
core flux density (Fig. ^ also seem to agree very well. The 
corresponding DCF analysis supports the correlation (Fig.|6}. 
There appears also a possible anti-correlation at about -250 
days, but this is very likely due to the broad peak around 2000.0 
overlapping with the dip at s/m2000.7 (Fig.0}. Calculation of 
the centroid of the highest 3 points of the DCF gives a short 
time delay of about 4 days for the spectral index, but at a bin 
size of 60 days the uncertainties are too large to reliably mea- 
sure such short delays. We thus performed cross-correlations 
between the optical R band light curve^ and our VLBI core 
light curves. The results are shown in Fig. The time sam- 
pling of the VLBI data is not as good as for the single-dish 
light curves. This might explain why the correlation is not so 
strong; however, all curves show a tentative correlation with a 
50 to 100 day time lag. The strongest correlation is found with 
the 43 GHz core light curve, which is also the best sampled one. 



3.3. Cross-correlation 

A cross-correlation of the single-dish light curves of BL Lac 

I ■ 1 

has already be en performed by Villata et al. ( 2004 a|) using a 

DCF analysis jEdelson & Kroliklil988i: iHufnagel & Breemanl 
I1992I) . We will continue to use this method here. The authors 
found that the higher frequency radio variations lead those 
at lower frequencies by several days to a few months, de- 
pending on the frequency separation. In particular, the varia- 
tions at the lowest frequency (5 GHz) lag the other frequen- 
cies by ~ 40 days (8 GHz), ~ 59 days (15 GHz), and ~ 78 
days (22 GHz) with uncertainties of 5 to 10 days. These val- 
ues were derived by calculating the centroid of the DCF near 
its peak. The corresponding uncertainties were estimated from 
variations of the time binning and the calculation of the cen- 
troid for cutoffs between 70% and 80% of the peak value. 
To complete the frequency coverage, we derive here also the 
DCF for the 5 GHz and 37 GHz frequency pai r and find a 
time lag of 105 + 15 days. To be consistent with lVillata et alJ 
ll2004a^ all DCFs in this work are performed for the same time 
range, between JD=2446500 (March 1986) and JD=2453000 
(December 2003), which excludes the highly pronounced ra- 
dio flare in 1980. Such prominent events could lead to spurious 
peaks in the DCF analysis, since a peak would appear at every 
time lag at which it overlaps with a flare in the second light 
curve. We note that the time range that we use is shorter for the 
VLBI data, which are only available between 1995 and 2004. 
All DCF calculations were performed several times with dif- 
ferent binning in time and over different time ranges. The plots 
shown here correspond to the versions that showed the most 
robust results. In most cases small changes did not aff'ect the 
results significantly and, where larger changes were observed, 
the se are discussed in th e text. 

IVillata et al.l (l2004al) found a modest correlation between 
the optical light curve and the high-frequency radio ones with 
a radio time lag of about 100 days. Moreover they showed that a 
better correlation is found when comparing the optical variabil- 
ity with the spectral variations in the radio bands, which high- 
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Fig. 6. DCF of the 43 GHz VLBI core flux density (1997-2002) 
and the radio spectral index between 8 GHz and 15 GHz show- 
ing a good correlation with a possible short delay of the spectral 
index. 

We obtain a better result from the cross-correlation of the 
optical data with the spectral indices. Figure |8] (top panel) 
shows the DCF for all four frequency pairs between 5 GHz and 
37 GHz (1986-2003). Although the DCF peaks are below 0.6, 
the clear trend for shorter time delays at higher frequency pairs 
is consistent with the radio cross-correlations, which show 
that the higher frequencies lead the lower ones ( Villata e t alJ 
2004a). The secondary peaks with negative time delays most 
likely come from the regular pattern of optical flares after 1 997, 
visible in Fig.|5] which appear at separations of a bit less than 
1 year (6 flares in 5 years). Auto-correlations of both the optical 
light curve and the 22/37 GHz spectral index weakly show this 
periodicity if we analyse only the time-range after 1997, but the 
periodicity vanishes outside of this 5-year period. Therefore, 

^ W e use optical flux densities corrected for the host galaxy contri- 
bution JViUata et all2002ll2004 j) 
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Fig. 7. DCF of the optical R band light curve with the VLBI 
core light curves at 15 GHz (1995-2003), 22 GHz and 43 GHz 
(1997-2002). 

we cannot rule out that the light curves of BL Lac sometimes 
show periodicities, but then it is only a short-term phenomenon. 
Periodicity analyses of BL Lac light curves performed by var- 
ious authors have found periods between 2.5 yr and 4 yr, about 
8 yr, and 1 5 yr to 2 yr for the radio (e.g.. Kellv et al. ,2003i 
fviUata et alJr2004ah . an d about 1 yr, 5yr, 7 to 8 yr, lOyr, and 
15 yrin the optical (e.g.. lHagen-Thorn et alJl 99711 Villata et alJ 
l2004ah . at different significance levels. It is evident that further 
investigations are needed before we can draw any conclusions 
regarding the presence of true, persistent periodicities. 

The formal statistical significance of the found time lags 
using the DCF analysis in this work is not easy to asses. In 
general the statistical significance of a correlation depends on 
the number o f points that co ntribute to the calculation of the 
time lag (e.g, IPetersonlEool . Since the number of points that 
contribute to the peak value of the DCF in our analysis always 
exceeds 100, all peaks higher than 0.25 have to be considered as 
being significant at a confidence level of about 98%. Therefore 
we carefully checked by eye how many flares with which am- 
plitude in the light curves contribute to the correlation and if 
there are any dominant flares which could strongly affect the 
DCF (see also Fig. Hence we might call the R - Qf22/37 
DCF with a peak value of 0.4 better than the correlation be- 
tween ag/is and the 43 GHz VLBI core light curve with a peak 
value 0.6, because the larger time interval and better sampling 
of R - a22/37 correlation makes it more reliable. A summary of 
the used parameters for all DCFs is given in Table|2 

Another method to test the robustness of the DCF peaks 
is to perform Monte Carlo simulations in which the curves 
are randomly modified within the measured uncertainties and 
only subsets of the data are correlated - a technique known 
as "flux redis tribution/rand om subset selectio n" (FR/RSS; 
IPeterson et ail 19981: see also lRaiteri et alJl2003l) . We used this 
method to test the peaks of the R — a cross-correlation (Fig.|8] 
(top)) through 1000 Monte Carlo simulations. In Fig. |8l (bot- 
tom) the fraction of occurrence of the 75% centroid values of 
all Monte Carlo realizations for the four DCFs are given. The 




Fig. 8. Top panel: DCF of the optical R band light curve with 
all spectral index pairs between 5 GHz and 37 GHz (1986- 
2003). Although the correlations reach a maximum peak value 
of only about 0.5, there is a clear trend to shorter time delays 
for higher frequencies. Bottom panel: Probability plot for the 
DCF peaks resulting from the correlation between optical data 
and radio spectral indices derived from 1000 Monte Carlo sim- 
ulations. 



histogram shows that the time delay and the uncertainty of the 
DCF between the R band and the radio spectral index decreases 
towards higher radio frequencies. 

To summarize: a fair correlation was found between 5 core 
and the radio spectral index with no measurable time delay, 
which suggests that the radio spectral index is a good indicator 
for the VLBI core variability. A tentative direct correlation is 
seen between S core and the optical light curve with radio time 
lags of less than 100 days. This, together with the correlation 
between the optical flux and the radio spectral index, provides 
evidence that the optical variability is indeed correlated with 
the radio emission of the VLBI core. 

In Table |3l we summarize the estimated time delays for the 
five radio frequencies during the time range of 1986 and 2003. 
Since the flattening of the spectrum appears when the higher 
frequency light curve rises earlier and/or faster than the lower 
frequency one, the variability of the spectral index should rep- 
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Table 2. Summary of cross-correlation parameters. Columns 
contain the frequency pairs (numbers denote radio frequencies 
in GHz and a "V" marks data obtained from VLBI images), 
time interval of correlation (MJD=JD-2400000), binning inter- 
val, DCF peak value, time-lag of the peak, time lag of the 75% 
centroid, and a reference to the corresponding figure. 
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resent the core variability at the higher frequency. Therefore, 
we attribute the time lag between the optical emission and the 
spectral index to the highest radio frequency. The delay of 
the 5 GHz variations, the lowest frequency light curve in the 
analysis, is derived by adding the time lag between the 5 GHz 
and 8 GHz radio light curves to the 7?-8GHz time lag. This 
procedure might lead to small shifts of the derived time lags, 
since the shape and peak of the spectral index variations are 
a combination of both frequencies, but tests with different fre- 
quency pairs that are more widely separated (e.g., 15/37 GHz or 
8/37 GHz) support this approach and yield shifts of less than 5 
days. This is well within the 1 cr errors derived from the Monte 
Carlo realizations. 

Table 3. Time lags, r, measured between the optical and the 
radio variations. The values are derived from the correlation of 
the optical emission with the radio spectral indices (see text for 
details). 
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Fig. 9. Optical light curve vs. radio spectral index (a8/i5)- The 
spectral index is shifted by -99 days. To make the comparison 
easier we represent the spectral index by (o-g/is -1-0.8)^ (top) and 
(ff8/i5 + 0.8)^ (bottom), so that all values are positive and the 
peaks are enhanced. The curve is smoothed with a three point 
average to reduce the noise. Top: All data. Bottom: Only the 
last 10 years with the best sampled data. 



Figure|9lillustrates the good agreement between the optical 
light curve and the radio spectral index (ag/is). In this figure 
the spectral index is shifted by -99 days according to the R- 
15 delay (TableO. Moreover, to make the comparison easier, 
we plot the spectral index by (ag/is -i- 0.8)"' (top) and (Qf8/i5 + 
0.8)'* (bottom), so that all values are positive and the peaks are 
enhanced. One can see that the curves agree fairly well during 
the last ~ 10 years (bottom panel). 

4. Discussion 

4.1. Optical-radio time delays 

The long observational history of BLLac led to several at- 
tempts to search for radio-optical correlations. Some of these 
found cor relations with radio t i me delays of several hun- 
dred days dPomphrev et alJll976t iHufnaeel & BreemarJI 19921: 



i Tornikoski eTaP Il994bt Iciements etaP 1 19951 'Hanski et al 
I2OO2'). Sometimes simultaneous events were reported ( Balonekl 
1982; Tornikoski et al. 1994a). A common result is that there 
are periods where radio and optical variations correlate better 
and others where the correlation seems to disappear. The corre- 
lations found here (Sect. lT3t are calculated for data from 1986 
to 2003, thus excluding the large outburst in 1980. Although the 
sampling is not homogeneous, it seems that the optical emis- 
sion underwent a period of suppressed activity from 1981 to 
early 1997 (Fig.|5Jl, during which flares were less frequent and 
not so pronounced. Such a period is visible neither in the radio 
light curves nor in the radio spectral index (Fig. 13. 

By constraining the DCF analysis on smaller time ranges 
we test how the time delays change with time. This was done 
exemplarily for the cross-correlation of the R band fluxes with 
a8/i5- The data between 1986 and 1989 are not so well sam- 
pled, but although the DCF results become less significant. 
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they confirm the previous findings. During the period before 
1986 a broad peak at about 300 days is found, which is in 
good agr eement with the results from the studies of data up 
to 1992 ("Pomphrev et al." 'l976' 'Hufnaeel & Breeman"l992; 
Tornikos ki et al. 1994b; Clements et al. 1995; Hanski et al. 
20021k during the period of weak optical variability (from 1986 
to early 1997) a small peak around zero time delay indicates 
some simultaneous events (Balonek 1982; Tornikoski et al. 
Il994al) . but in general the DCF during this time shows no peak 
above 0.25. Finally the most recent data (1997 to 2003), char- 
acterized by the high optical activity, are d ominated by a p eak 
at around 100 days. The fact that both Vill ata et al ] |l2004ah and 
our new analysis detected this 100-day time delay for data ob- 
tained from 1986 until now, demonstrates that the most recent 
correlation seems to dominate over that found in the previous 
period of weak variability. However, the time range used here 
is not dominated by a major outburst, but rather consists of sev- 
eral flares in the optical (6-7 major events after 1997 and about 
the same number of minor flares before 1997), and the vari- 
ability of the spectral index displays nearly the same amplitude 
throughout the whole period. Therefore, the optical-radio de- 
lays found in Sect. 13.31 seem to represent the result of a true 
correlation in the variability behaviour of BLLac over at least 
the last 10 years of our dataset. 



4.2. Frequency dependence of the delay 

Frequency dependent time delays in radio flares are commonly 
related to a combination of optical depth and travel time along 
the jet. Evolving flares first appear in the inner jet region at 
the highest frequencies and, as the perturbations (VLBI com- 
ponents) travel down the jet, the flares become visible at suc- 
cessively lower frequencies where the optical depth of syn- 
chrotron self-absorption decreases to Tj ~ 1 at that frequency. 
For a given magnetic field strength (B(r)) and electron energy 
distribution sca le factor (N{r)), the corresponding Ts is (e.g., 
lT,obanovll998h: 



C2(a)Ni 



eBi 
2nm, 



^{em+;i-l)y6+l 



(1) 



Here nig is the electron mass, e is the electron charge, 
e - 3/2 - a, a is the optically thin synchrotron spectral index, 
and m and n are the power-law exponents of the radial distance 
dependence of the magnetic field and N{r), respectively (e.g., 
lLobanovll99"8l) . The observed jet opening angle is 0o - (p esc 6. 
The factor €2(0:) is described in Blumenthal & Gould ( 1970), 
and C2(-0.5) - 8.4 • 10'° in cgs units. Setting Ts to unity gives 
us an estimate of the distance from the convergence point of 
the jet to the core as observed at frequency v. 

Hence, we should be able to measure a frequency- 
dependent core position in VLBI images: for a conical jet ge- 
ometry the shift is giv en by r oc v~^^^% where ky - [(3 - 2a)m + 
2n - 2]/(5 - 2a) (e.g.. iLobanovl 1 99^ . For the case of equipar- 
tition between the magnetic field and electron energy densi- 
ties, the simplest choice of m = 1 and n - 2 (Konisl 1981) 
leads to = 1 independent of the spectral index dLobanovl 
Il998l) . Larger values are reached in regions with steeper pres- 



sure gradients than if, e.g., the electrons suffer adiabatic en- 
ergy losses without reacceleration (Marscher 1980). Thus, tak- 
ing measurements of the core shift between several frequen- 
cies allows to estimate A:, along the jet (e.g.. lLobanovlfl995 
Kadleret al. 2004). 

If the radio time lags are also due to opacity effects and 
the flare travels at a constant speed, then they should also be 
proportional to v . Recent kinematic studies on high fre- 
quency VLBA data of BL Lac have consistently derived jet 
speeds with a bulk Lorentz factor of about 7 and an angle 
to the line of sight of ~ 12°, with only minor deviations 
jPenn et alJ boOOf iRevnolds et alJ l2003t IStirling et all l2003t 



'kellermann et al."2004'; 'jorstad et al."2005V We therefore as- 
sume that the Lorentz factor and angle to the line of sight are 
more or less constant, and that the time delays are proportional 
to the distance between the two emitting regions. E.g., a time 
delay, f', of 10 days would then correspond to a distance, r, of 

about 0.06 pc (r - t(3c, with /3 - ■yj'^ - ~t and t - yf ). 

In Fig.[TO|we plot the measured radio lags with respect to 
both the 37 GHz data and to the optical data vs. frequency. A 
power-law fit represents the data reasonably well in both cases 
(reduced < 1) and there it seems that the assumption of a 
constant jet speed holds at least in the radio regime. It is im- 
portant to note that the value of strongly depends on the ab- 
solute time lag, i.e. on the distance to the base of the jet, which 
might be not far from the central engine. Although the optical 
emission could be emitted close to the foot point of the jet, we 
should consider that this region is still at some considerable dis- 
tance from the central engine. In this case A:, would be larger. If, 
on the other hand, the jet speed between the radio and the op- 
tical emitting regions is not constant, the simple assumption of 
T oc v^'/'^f is invalid. However, it seems that at least in the radio 
regime the speed is constant and we can infer that k^ is likely 
larger than 1.1. This implies that the gradients in jet parameters 
are steeper than in the case of equipartition with m - I. Our re- 
sult is, on the other hand, consistent with a freely expanding jet 
with adiabatic en ergy losses, in which case n - (2/3)(3 - 2a) 
( lMarscherll980t ) and k, = (15 - 14a)/[3(5 - 2a)] for m = 1, or 
~ 1.2 for a ^ -0.5. 



4.3. Connection with jet components 

Since outbursts (especially at mm wavelengths) are often re- 
lated to the ejection of superluminal jet c omponents (e.g, 
lOtterbein et'ai]ll998t ISavolainen et al.ll2002l) . we have com- 
pared our radio and optical light curves and the spectral index, 
ag/is, with the ejection dates of jet components reported in the 
literature. Figure^JiUustrates this comparison. There is no ob- 
vious connection among the events, but some new components 
are led by an optical flare or followed by an increase in the ra- 
dio. We also tested whether we could find any evidence for a 
connection between the ejections and flares by calculating the 
variations in the flux density at different frequencies and the 
spectral indices for times shortly before and after the ejection, 
but no clear trend could be found. There is always a significant 
fraction of the components that does not show the expected be- 
haviour. However, this kind of analysis might be affected by the 
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Fig. 11. Top: Optical (grey squares) and radio (black circles) light curves in comparison with the ejection dates of new super- 
luminal jet components, which are highlighted as grey vertical bars. The width of each bar corresponds to the uncertainty of 
the ejection date. Bottom: Ejection dates vs. radio spectral index (b lack circles and gr e y trend-line) and the 43 GHz core flux 
density (black squares and line). The ejection dates are taken from* Denn et al.l (l2000l) . ]s'tirling et al.l (l2003l) (S7 to S13), and 
iMarscher et al. (.2004.) (Ml and M2). Only some of the new components seem to be connected with events in the light curves. 



^Rband ' 
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l:,=l.S3 + 0.13 
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Fig. 10. Measured time lags vs. frequency. The circles and 
dashed line represent the measured time lags and the corre- 
sponding fit with respect to the optical emission, while the 
squares and dotted line denote time lags and fit with respect 
to the 37 GHz light curve. 



relatively large errors in the dates of the earlier ejections (S7 
to SIO) and the fact that some components appear very close 



together (S9/S10 and S12/S13). Furthermore, the light curves 
are generally composites of multiple events, some of which are 
declining in flux while others are rising. This can render indi- 
vidual flares unnoticed in the total-flux light curves. 

The best agreement can be seen for the two newest com- 
ponents Ml and M2. They appear accompanied by an optical 
flare and are followed by a flare and flattening of the spectrum 
at radio frequencies. Since previous works have also reported a 
good agreement between the appearance of flares and the ejec- 
tion of com ponents from SI up to S7 (e.g., iMutel et 3111199(1 
|Aner"l999), it might be that the dates of S8 to S13 are not so 
well determined. Although S7, Sll, and S12 can be identified 
with minor radio flares and a flattening of the spectrum, there 
are several large flares in the optical and radio light curves at 
1997.6, 1998.4 and 1999.5 that seem to be unrelated to any new 
component. It will be interesting to see if the ongoing densely 
sampled VLBI monitoring will find corresponding components 
to the pronounced optical flares around 2001.5 and 2002.1. 

As mentioned earlier, one should also expect to see a sign 
of the travelling jet components in the different parts of the 
VLBI light curves. After a flare in the VLBI core, a new com- 
ponent should brighten the jet after some time at every mea- 
sured distance while it passes by. Unfortunately, the jet flux 
density is either too low or the time sampling is insufficient 
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to reliably detect such events. In addition, the ejection of sev- 
eral components or the appearance of many flares in a short 
time can confuse the picture. A good example of the evolu- 
tion of a new jet component is Ml (Fig. II 1> ejected in 1999.95 
(JD=2451527). The corresponding flare is a nicely separated 
event and from Fig.|3]one can see that, some time after the peak 
in the core, the inner jet (r x 0.4 mas) brightens at 43 GHz and 
22 GHz, while later it appears at 2 mas and finally at 3 mas sep- 
aration from the core. Since not all frequencies were observed 
simultaneously and in order to see better the evolution, we have 
converted the flux densities at 15 GHz and 43 GHz of the core 
and jet to 22 GHz according to their measured spectral index 
and show these together in Fig. El The plotted data points are 
three-point averages to smooth the residual differences, and the 
error bars correspond to their standard deviations. Parabolic 
fits to the data illustrate the appearance of the component at 
the diff'erent jet positions. The calculated speed of the compo- 
nent from the peaks in the Ught curves is about (2-4) mas/yr or 
jSapp ~ (12-19)c. This is somewhat higher than the speeds of 
5 c to 9 c from kinematic studies reported thus fa r (Denn et al. 
2000; St irling et alJ2003UKellermann et alJ20o4ijorstad et alJ 
2005), but since we have combined all frequencies in order to 
be able to fit the peak, our values might be more uncertain. 
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Fig. 12. VLBI light curves of the radio flare in early 2000. To 
better follow the flare evolution in both time and space, all 
fluxes are converted to 22 GHz according to their spectral in- 
dex and combined in a three point average. One can follow 
the propagation of the core outburst along the jet. The ejec- 
tion dates of new VLBI components are marked by arrows and 
correspond to MI in 1999.95 and M2 in 2000.81 reported by 
lMarscheretaD ( l2003l) . 

Another jet comp onent ejection (M2) is reported for 
2000.81 (JD^245 1 841 : iMarscher et"al]l2003l) and the related 
radio flare is also visible in Fig.^] There is a sharp core peak 
at JD=2451920 and some time later a peak in the r « 0.4 mas 
light curve. Unfortunately, the dense VLBI monitoring ends in 
early 2001 and one cannot follow the flare further. However, 
both ejections are accompanied or led by a bright optical state, 
which suggests that the correlated optical-radio variability ob- 
served here is also related to the ejection of jet components. 



During 1999 and 2000 BLLac was the target of several 
j^Xr^ ( Mar scher et al. 2004) an d BeppoSAX iRavasio et ali 
I2OO2I r2003t iBottcher etalJl2003l) observations. Interestingly, 
the source was caught in very different states. The RXTE light 
curve covers nearly the whole period of 1999 to 2001 with 
dense sampling. In particular, during the ejection in 2000.81 
strong variability is observed together with a softening of the 
X-ray spectrum (Marscheret al. 2004). This is also confirmed 
by the BeppoSAX observ ations from October 3 1 to November 
2, 2000 (JD~2451850; iRavasio et alJ l2003[ iBottcher et all 
2003), which show an unusually strong soft spectral compo- 
nent between 0.2 keV and about lOkeV. Unfortunately, the 
flare/ejection in 1999.95 was not covered by RXTE, while 
the BeppoSAX observations of December 5-6 (JD~2451518) 
obtained several days before the extrapolated ejection date 
showed no sign of a soft component, but a hard spectrum simi- 
lar to that observed in July 2000. Another BeppoSAX observa- 
tion on June 5-7, 1999 (JD~2451336) shows aga in a small soft 
component between 0.2 keV and about 3 keV ^Rav asio et alJ 
2002), which seems also visible as^ s mall enhancement in the 
RXTE light curve and spectrum dMarscher et alJl2004l) . This 
small X-ray flare is also followed by optical and radio flares, 
but no new jet component has been reported. However, this 
might be due to confusion with the two earlier ejections of com- 
ponents S12 and S13, which were separated by only 0.3 yr. 

Although further coordinated observations are needed to 
complete the picture, the very well sampled period from 1999 
to 2001 already reveals some very interesting similarities in the 
variability behaviour from X-ray to radio wavelengths, which 
are also visible as changes in the source VLBI structure. The 
missing soft-spectrum component in the BeppoSAX observa- 
tions of 1999.92 might blear the claim of a correlation, but 
given that the typical uncertainty in the ejection dates is about 
0.2 yr and the soft X-ray tail is potentially a short-lived phe- 
nomenon, more observations yielding better statistics are nec- 
essary. 

4.4. A possible jet scenario 

Due to the combination of the structural analysis using VLBI 
light curves and the single-dish observations, we are con- 
fident that the emission from the VLBI core, presumably 
close to the base of the jet, is responsible for the corre- 
lated variability observed at optical and radi o wavelengths, 
and possibly also up to X-ray frequencies (*Mars cher et alJ 
|2004). The mechanisms that can explain blazar variability 
are still the subject of discussion. P ossibilities suggested in 
the literature include shocks in jets (" Marscher & Geajll985t 
pVller et al. 1985; Marscher 1996), and flares introduced by 
the lighthouse effect from changes in the direction of forward 
beaming caused, e.g., by helical trajectories of plasma ele- 
ments (Camenzind & Krockenberget 1992), by a p recessing bi- 
nary b lack-hole system t .Begelman et al . 1980; Sillanpaa et al.) 
11988"), or by flie rotation of a he lical jet tVillata et al...l99&: 
Villata & Raiteril [l999; Ostorer o et all l2004 . Furfliermore. 
variability cause d by colliding relativistic p lasma shells has 
been considered dSpada et alJ200ltlGuetta et al.i.2004) . 
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Regardless of the variability origin, the VLBI observations 
of BLLac indicate a bent path of the jet. The trajectories of 
several jet components have been well modelled and predicte d 
by helical jet models JPenn et alJbOOOt IStirling et al] l2003'). 
A bent jet path could also explain the region of low emission 
around 1 mas from the core noticed in Sect. 12.21 (Fig. l2t. One 
can see strong variability in both the core and the jet at core 
separations larger than 1 mas, but never between 0.7 mas and 
1 mas. Since the spectral index does not show any evidence 
for absorption, a natural explanation would be reduced Doppler 
boosting either from a misalignment between the jet direction 
and our line of sight or from a change in the jet speed. 

A possible scenario that could account for the observed 
variability behaviour, especially the changing variability pat- 
tern between optical and radio emission on timescales of sev- 
eral years, would be a helical jet structure that is precessing 
or rotating. Since the optical emission is most likely produced 
on small scales, a small misalignment could beam the emission 
less favourably, which could explain the low optical state with 
suppressed variability from 1981 to early 1997 (see Sect. 14. 1> . 
On the contrary, the radio emission, coming from all along the 
jet, would be always strongly beamed in some part of the heli- 
cal path and therefore would not suffer t he less intens e beaming 
of the optical component (e.g., lRiegerll2004 (20051 and refer- 
ences therein). The only beaming-suppressed region at radio 
frequencies would be the emission gap in the inner jet. In this 
case the gap might move or disappear with time and could be 
used by future VLBI studies to prove or disprove the rotating 
jet scenario. 

5. Summary 

We have presented an analysis of a VLBI data set that includes 
108 images at 15 GHz, 22 GHz, and 43 GHz obtained between 
1995 and 2003. The aim has been to disentangle the different 
components contributing to the single-dish radio light curves. 
Inspection of the VLBI light curves corresponding to different 
regions of the jet reveals that the radio single-dish light curves 
indeed display the sum of the emission from the core and from 
some prominent jet features. In some cases these jet features 
become as bright as, or even brighter than, the core itself. This 
leads to a complication of the cross-correlation between the ra- 
dio light curves and other wavebands where the emission might 
come from the core region alone. 

We find that the spectral variability that is present in the 
single-dish light curves can trace the variability of the VLBI 
core, and therefore enables a study of the long-term variability 
of the radio core even when VLBI data are not available. Using 
this result, we find a fair correlation between the variations of 
the radio spectral index and the /?-band optical light curve with 
radio delays of about 50 to 180 days depending on the fre- 
quency separation. (Here we attribute the delay measured from 
the spectral index between two frequencies to the higher fre- 
quency. This seems justified as the flattening and also the peak 
of the spectral index curve is mainly due to the earlier and faster 
rise of the flare at the higher frequency.) Owing to the shorter 
time range and the sparser sampling, the radio-optical correla- 
tion is weaker when considering the VLBI core light curves at 



22 GHz and 43 GHz directly, but the delays found are compa- 
rable to those of the spectral indices. The resulting time lags 
are also c onsistent with those f rom the cross-correlations per- 
formed by'Villata et al.' ('2004aV 

Assuming that the radio time lags are due to an optical 
depth effect of synchrotron self-absorption and that flares prop- 
agate along the jet at a constant speed, we obtain a power-law 
dependence of the time lag on frequency that suggests that 
the jet does not maintain equipartition of magnetic and parti- 
cle energy densities as it expands. The radio-optical correlation 
is most prominent before 1981 and after 1997. During the 16 
years in between, we observe only moderate variability in the 
optical and the correlation with the radio almost disappears. 
Stirling et al. (2001) have suggested a precessing jet nozzle for 
BL Lac to explain their observations, which may imply that the 
optical variability is also orientation dependent. This in turn, 
depending on the geometry, could also affect the correlation 
with the radio regime. This would lead to a natural explana- 
tion of the changing appearance of correlated variability: since 
the optical emission likely originates on small scales, a minor 
misalignment could beam the emission less favourable. On the 
other hand, the beamed radio emission could still repeatedly 
approach our line of sight along its helical path farther down 
the jet, so the precession or rotation might modify the ampli- 
tude of the variability pattern without leading to a vanishing 
flux as in the optical bands (e.g., Rieger 2004, 2005). A bent 
jet path could also explain the low emission region in the jet 
around 1 mas from the core. A movement or the disappearance 
of the gap in future VLBI studies, could therewith be used to 
prove or disprove a rotation of the bent jet path. 
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Table .1. Observing log. Given are the observing epoch, frequency, beam size and orientation, total flux in the image (Stot), peak flux density 
(Smax). rms noise level (Snns). and the references where the data was first published. 



Epoch 


V 


Beam 




Stot 


Smax 


Srms 


Ref. 




[GHz] 


[mas X mas] 


[°] 


[Jy] 


[Jy] 


[mJy] 




1995.27 


15 


0.50x0.88 


2.9 


3.30 


1.79 


2.4 


(1) 


1995.96 


15 


0.53 X 0.97 


0.9 


5.01 


2.04 


2.5 


(1) 


1996.38 


15 


0.51 X 0.90 


-4.1 


5.67 


2.68 


2.4 


(1) 


1996.82 


15 


0.46 X 0.86 


-2.6 


3.84 


1.77 


3.2 


(1) 


1997.19 


15 


0.58x1.00 


7.1 


3.21 


1.49 


1.5 


(1) 


1997.66 


15 


0.56x0.91 


12.4 


2.97 


1.70 


1.5 


(1) 


1998.18 


15 


0.50 X 1.00 


8.2 


2.68 


1.06 


3.8 


(1) 


1998.73 


15 


0.57 X 0.87 


-16.0 


4.07 


2.48 


0.7 


(5) 


1998.97 


15 


0.59x0.71 


-3.0 


2.98 


1.54 


0.7 


(5) 


1999.16 


15 


0.63 X 0.78 


6.5 


3.26 


1.88 


0.6 


(5) 


1999.40 


15 


0.60 X 0.81 


6.2 


2.38 


1.30 


0.7 


(5) 


1999.70 


15 


0.55 X 0.77 


-14.1 


2.75 


1.71 


0.6 


(5) 


1999.85 


15 


0.59 X 1.08 


-3.3 


3.31 


2.34 


1.3 


(1) 


2000.01 


15 


0.56 X 0.79 


-16.1 


4.23 


3.32 


0.7 


(5) 


2000.31 


15 


0.55 X 0.74 


-4.7 


4.18 


3.16 


0.6 


(5) 


2000.49 


15 


0.56 X 1.00 


16.7 


3.23 


1.96 


0.5 


(4) 


2000.84 


15 


0.57 X 0.97 


27.1 


2.14 


1.32 


0.4 


(4) 


2000.99 


15 


0.54 X 0.94 


-16.3 


2.59 


1.82 


1.1 


(1) 


2001.00 


15 


0.47 X 1.46 


-6.0 


2.83 


2.08 


0.8 


(4) 


2001.09 


15 


0.49x1.27 


-10.9 


2.46 


1.79 


1.1 


(4) 


2001.17 


15 


0.52 X 1.37 


-1.8 


1.92 


1.35 


0.9 


(4) 


2001.22 


15 


0.57 X 1.02 


33.8 


1.78 


1.14 


0.3 


(4) 


2001.25 


15 


0.54x1.44 


-1.1 


1.96 


1.26 


0.6 


(4) 


2001.97 


15 


0.48 X 0.92 


1.6 


2.34 


1.44 


0.7 


(1) 


2002.03 


15 


0.32 X 0.64 


-33.2 


2.62 


1.40 


0.4 


(6) 


2002.05 


15 


0.85x1.15 


29.3 


2.75 


1.79 


0.6 


(5) 


2002.45 


15 


0.52 X 0.79 


-18.1 


2.18 


1.24 


0.7 


(2) 


2002.51 


15 


0.24 X 0.55 


-20.4 


2.52 


1.38 


0.5 


(6) 


2003.10 


15 


0.51 X 0.76 


-0.1 


1.75 


0.98 


0.5 


(2) 


2003.24 


15 


0.29 X 0.46 


-35.8 


1.73 


0.80 


0.3 


(6) 


2003.65 


15 


0.48x0.71 


0.8 


3.11 


2.40 


0.7 


(2) 


2004.44 


15 


0.53 x 0.73 


-9.3 


2.73 


2.01 


0.7 


(2) 


1997.94 


22 


0.39x0.84 


-14.1 


2.79 


1.54 


1.2 


(4) 


1998.03 


22 


0.45 X 0.87 


-8.0 


3.06 


1.51 


1.1 


(4) 


1998.10 


22 


0.34 X 0.85 


-11.3 


2.81 


1.21 


1.3 


(4) 


1998.18 


22 


0.36 X 1.00 


1.9 


2.90 


1.42 


1.6 


(4) 


1998.27 


22 


0.35 X 0.85 


-8.2 


6.60 


3.36 


2.2 


(4) 


1998.35 


22 


0.36x0.83 


-7.4 


3.17 


1.58 


1.5 


(4) 


1998.44 


22 


0.35 X 0.87 


-8.6 


2.79 


1.48 


1.4 


(4) 


1998.52 


22 


0.33 X 0.97 


-13.6 


3.22 


1.82 


1.5 


(4) 


1998.61 


22 


0.36x0.87 


-10.4 


3.97 


2.71 


2.1 


(4) 


1998.73 


22 


0.42 X 0.64 


-17.8 


3.25 


2.13 


1.0 


(5) 


1998.82 


22 


0.35 X 0.93 


-21.4 


3.30 


1.88 


1.4 


(4) 


1998.92 


22 


0.36x0.88 


-14.8 


2.72 


1.44 


1.5 


(4) 


1998.97 


22 


0.48 X 0.57 


-5.9 


2.52 


1.38 


0.6 


(5) 


1999.03 


22 


0.38 X 1.02 


-13.2 


3.72 


2.07 


1.1 


(4) 


1999.11 


22 


0.36x0.82 


-8.3 


2.82 


1.79 


1.6 


(4) 


1999.16 


22 


0.44 X 0.56 


1.6 


2.83 


1.72 


0.5 


(5) 


1999.41 


22 


0.41 X 0.58 


-5.5 


1.93 


1.06 


0.6 


(5) 


1999.70 


22 


0.39 X 0.64 


-26.4 


2.07 


1.37 


0.7 


(5) 


2000.01 


22 


0.43 X 0.59 


-15.7 


4.32 


3.43 


0.5 


(5) 
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Table . 1 - continued from previous column 



Epoch 


V 


Beam 




Slot 


Smax 


Srms 


Ref. 




[GHz] 


[mas X mas] 


n 


[Jyl 


[Jyl 


[mJy] 




2000.31 


22 


0.43 X 0.56 


-8.1 


3.76 


2.76 


0.5 


(5) 


2000.49 


22 


0.40 X 0.74 


12.8 


3.20 


2.02 


1.3 


(4) 


2000.84 


22 


0.45 X 0.69 


31.4 


2.20 


1.43 


0.8 


(4) 


2001.00 


22 


0.37 X 0.96 


-4.0 


2.65 


2.07 


1.3 


(4) 


2001.09 


22 


0.38x0.84 


-13.3 


2.38 


1.65 


1.5 


(4) 


2001.17 


22 


0.40 X 0.94 


-1.9 


1.69 


1.20 


1.1 


(4) 


2001.22 


22 


0.46 X 0.78 


33.7 


1.98 


1.34 


0.7 


(4) 


2001.25 


22 


0.39x0.98 


-4.8 


1.86 


1.21 


0.8 


(4) 


2001.59 


22 


0.37x0.63 


-11.3 


1.88 


1.34 


0.7 


(5) 


2002.05 


22 


0.65 X 0.92 


24.0 


2.28 


1.61 


0.9 


(5) 


1997.86 


43 


0.17x0.43 


-10.5 


3.67 


2.74 


2.1 


(4) 


1997.94 


43 


0.18 X 0.44 


-10.2 


1.75 


0.98 


1.8 


(4) 


1998.03 


43 


0.23 X 0.46 


-9.1 


2.02 


1.04 


1.1 


(4) 


1998.10 


43 


0.17x0.43 


-12.2 


1.93 


0.97 


1.6 


(4) 


1998.18 


43 


0.17x0.53 


-2.7 


2.30 


1.33 


1.8 


(4) 


1998.23 


43 


0.17x0.33 


9.9 


3.49 


1.96 


1.0 


(3) 


1998.27 


43 


0.18 X 0.44 


-10.2 


3.56 


2.17 


1.9 


(4) 


1998.35 


43 


0.16x0.44 


-10.5 


2.35 


1.30 


2.4 


(4) 


1998.41 


43 


0.16x0.32 


-14.6 


2.70 


1.35 


0.6 


(3) 


1998.44 


43 


0.18x0.44 


-9.4 


2.01 


1.13 


2.1 


(4) 


1998.52 


43 


0.17x0.57 


-15.6 


3.21 


1.96 


1.8 


(4) 


1998.58 


43 


0.16x0.28 


1.1 


4.28 


2.98 


0.8 


(3) 


1998.61 


43 


0.19x0.43 


-10.6 


3.85 


2.80 


1.8 


(4) 


1998.70 


43 


0.17x0.56 


-17.3 


3.88 


2.65 


1.4 


(4) 


1998.73 


43 


0.20 X 0.32 


-25.9 


3.07 


1.90 


0.9 


(5) 


1998.76 


43 


0.19x0.30 


-11.5 


2.50 


1.44 


0.7 


(3) 


1998.82 


43 


0.17x0.47 


-19.3 


2.33 


1.34 


1.7 


(4) 


1998.92 


43 


0.21 X 0.49 


-23.6 


1.99 


1.18 


1.4 


(4) 


1998.94 


43 


0.19x0.30 


-13.4 


2.00 


1.14 


0.7 


(3) 


1998.97 


43 


0.21 xO.26 


-11.3 


1.86 


0.99 


0.7 


(5) 


1999.03 


43 


0.17x0.60 


-11.8 


2.96 


1.71 


1.6 


(4) 


1999.11 


43 


0.17 X 0.42 


-9.7 


2.05 


1.47 


1.7 


(4) 


1999.12 


43 


0.19x0.32 


-13.9 


2.71 


1.64 


0.6 


(3) 


1999.16 


43 


0.21 X 0.27 


-1.1 


2.31 


1.38 


0.6 


(5) 


1999.32 


43 


0.18x0.26 


7.8 


2.10 


1.19 


1.2 


(3) 


1999.41 


43 


0.23 X 0.39 


-18.7 


1.50 


0.84 


0.7 


(5) 


1999.70 


43 


0.19x0.35 


-30.0 


2.12 


1.50 


0.7 


(5) 


1999.76 


43 


0.20 X 0.38 


-16.8 


2.73 


2.00 


1.2 


(3) 


1999.93 


43 


0.20x0.38 


-15.9 


3.68 


2.76 


1.0 


(3) 


2000.01 


43 


0.23 X 0.31 


-9.7 


3.34 


2.58 


0.9 


(5) 


2000.06 


43 


0.19 X 0.26 


7.2 


4.21 


3.00 


1.7 


(3) 


2000.26 


43 


0.18 X 0.28 


1.5 


4.17 


2.45 


1.0 


(3) 


2000.31 


43 


0.19x0.26 


2.0 


2.95 


1.80 


0.9 


(5) 


2000.49 


43 


0.21 X 0.37 


20.7 


2.39 


1.37 


1.6 


(4) 


2000.54 


43 


0.18 X 0.30 


-6.3 


2.03 


0.97 


0.8 


(3) 


2000.75 


43 


0.19x0.32 


-8.8 


1.47 


0.74 


0.8 


(3) 


2000.84 


43 


0.23 x 0.35 


32.5 


1.88 


1.23 


1.7 


(4) 


2000.94 


43 


0.23 X 0.33 


28.3 


1.73 


1.19 


1.1 


(3) 


2001.00 


43 


0.18x0.49 


-4.6 


2.83 


1.98 


1.2 


(4) 


2001.09 


43 


0.18 X 0.44 


-12.6 


1.94 


1.18 


1.7 


(4) 


2001.17 


43 


0.17x0.49 


-4.4 


1.46 


0.92 


1.5 


(4) 


2001.22 


43 


0.21 X 0.36 


33.4 


1.61 


0.99 


1.2 


(4) 


2001.25 


43 


0.19x0.48 


-3.7 


1.63 


1.02 


1.5 


(4) 


2001.28 


43 


0.17x0.29 


-2.6 


1.44 


0.83 


0.6 


(3) 
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Table . 1 - continued from previous column 



Epoch 


V 


Beam 




Stot 


Smax 




Ref. 




[GHz] 


[mas X mas] 


[°] 


[Jy] 


[Jy] 


[mJy] 




2001.59 


43 


0.26 X 0.46 


13.1 


1.99 


1.48 


1.0 


(5) 


2002.05 


43 


0.33 X 0.57 


33.3 


1.73 


1.27 


1.4 


(5) 


2003.82 


43 


0.20x0.38 


-22.1 


1.88 


1.27 


0.7 


(5) 



Note: ( 1 ) VLBA 2 cm Survey JKellermann et all 1 1 998t IZensus et al J Eooi iKellermann etal] l2004 : (2) MOJAVE 
jLister & HoniarJl2005l): (31 1.Torstad et al.1 (l200l[ l2005l) . LStirling et alJ hm^ : (4) Gomez et al. & Agudo et al. priv. comm.; 
(5l lDenn et alJll2000l) . lMutel & Dennlll2005l) : (61 this work. 
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